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Ground State Properties of One Dimensional S = 1/2 Heisenberg 
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The one dimensional S = 1/2 Heisenberg model with dimerization and quadrumerization is 
studied by means of the numerical exact diagonalization of finite size systems. Using the phe- 
nomenological renormalization group and finite size scaling law, the ground state phase diagram 
is obtained in the isotropic case, ft exhibits a variety of the ground states which contains the 
5 = 1 Haldane state, 5 = 1 dimer state and 5 = 1/2 dimer state as limiting cases. The gap 
exponent v is also calculated which coincides with the value for the dimerization transition of 
the isotropic Heisenberg chain. In the XY limit, the phase diagram is obtained analytically and 
the comparison is made with the isotropic case. 
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§1. Introduction 

Recently, the S = 1/2 antiferromagnetic Heisenberg 
chains (AFHC) with modulated spatial structures have 
attracted a great deal of attention. Although the uni- 
form 5 = 1/2 AFHC can be solved exactly using the 
Bethe Ansatz^ method and is known to have a gapless 
ground state,EP this state is unstable against the dimer- 
ization leadintto the spin-Peierls state which has the 
spin gap.DBclEP Other kinds of spatial modulation of 
the exchange coupling arising from the underlying lat- 
tice structures can also induce various kinds of spin gap 
phases. 

On the other hand, no spin-Peierls instability is ex- 
pected in tha.S' = 1 AFHC due to the presence of the 
Haldane gap.B' Instead, the Haldane-dimer phase ti; 
sition takes place at finite strength of dimerization! 
One of the authors (KH) has pointed out the connec- 
tion between the dimer state of the 5 = 1/2 dimerized 
Heisenberg chain and the Haldane state of the 5 = 1 
AFHC.lP In this context, the dimerization in the 5 = 1 
chain corresponds to the quadrumerization in the dimer- 
ized 5 = 1/2 chain. From this point of view, in the 
present work, we aim to get deeper insight into the na- 
ture of the Haldane-dimer transition in the 5=1 AFHC 
from the investigation of the ground state properties of 
the 5 = 1/2 Heisenberg chain with dimerization and 
quadrumerization. 

On the other hand, in the neighbourhood of the 5 = 
1/2 uniform Heisenberg point, both dimerization and 
quadrumerization are expected to produce the energy 
gap. It is possible, however, that the coexistence of these 
two periodicities does not always enhance the gap addi- 
tively but may rather reduce the gap due to the competi- 
tion between them leading to the gapless state in spite of 
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the spatial nonuniformity. We also intend to investigate 
such competition of two periodicities in this model. 

This paper is organized as follows. In the next sec- 
tion, the model Hamiltonian is defined. The numerical 
results are presented in §3 for the isotropic case. From 
the phenomenological renormalization group and finite 
size scaling analysis of the numerically calculated energy 
gaps, the ground state phase diagram is obtained in §3. 
The critical exponent v of the energy gap is also esti- 
mated taking the logarithmic corrections into account. 
The numerical results are compared with the perturba- 
tion theory. In §4, the ground state phase diagram is 
obtained analytically for the XY case using the Jordan- 
Wigner transformation. The comparison is made with 
the isotropic case. The last section is devoted to sum- 
mary and discussion. 

§2. Model Hamiltonian 

The Hamiltonian of the one dimensional dimerized and 
quadrumerized 5 = 1/2 Heisenberg chain is given by 

2N 

= ^ [j('5'2i-l'5'2i + ^21-1^21 + ^^2l-1^2l) 
1 = 1 

2N-1 



^ [(l + (-l)'-i5)(5f,5; 



21+1 



qV qV 
'^2l'^2l+l 



1=1 



^^2l^2l + l] 



(2.1) 



where 1 — i{—oo < j < oo) and d{—l < 6 < 1) represent 
the degree of dimerization and quadrumerization, respec- 
tively. The open boundary condition is assumed. The 
anisotropy parameter is denoted by A. In the present 
work, we concentrate ourselves on the cases A = 1 
(isotropic case) and {XY case). We also take 6 > 
without loss of generality. 

Let us discuss the ground state of this model in some 
limiting cases for A = 1. First, our model tends to the 
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S = 1 AFHC for large negative j. Therefore the ground 
state is the VBS-hke Haldane phase or the S — 1 dimer 
phase according as ^ < (5c or 5 > 5c where Sc is the crit- 
ical vabjfi which depends on j. For j —oo, Sc tend 
to 0.25lI'I3) which is the value for the S* = 1 AFHC. In 
terms of the spin-1/2 language, these phases can be also 
described as follows. For small S (Haldane phase), the 
spin pairs connected by the 1 + 5-bonds and 1 — 5-bonds 
form local singlet pairs as schematically shown in Fig. 
1(a). For large S {S — I dimer phase), the spins con- 
nected by 1 -I- (5-bonds {84,1+2 and 84,1+3) form singlet 
pairs strongly. Mediated by the fluctuation within these 
singlet pairs, the effective antiferromagnetic coupling is 
induced between the spins S41+1 and S41+4, which leads 
to the 4-spin local singlets as shown in Fig. 1(b). It 
should be noted that the above picture in spin-1/2 lan- 
guage remain valid even for positive j as far as j < 1. 

Secondly, we discuss the neighbourhood of S — and 
J ~ 1. For (5 = 0, the singlet pairs reside on the j-bonds 
or the 1 ± (5 bonds according as j > 1 or < 1 correspond- 
ing to the 5' = 1/2 dimerized state with oppsite parity. 
For J < 1, the dimer configuration is the same as in Fig. 
1(a) while it is shown in Fig. 1(c) for j > 1. 

§3. Numerical Results for A = 1 

The Hamiltonian is numerically diagonalized to cal- 
culate the energy gap G{N,j,S) in the open chain for 
AN ~ 12, 16, 20 and 24 using the Lanczos algorithm. The 
open chain is used so that the critical point can be eas- 
ily discerned. For concreteness, let us take the leftmost 
and rightmost bonds as j-bonds. According to the spin 
structure discussed is §2, for small values of S, there re- 
main two S* = 1/2 residual spins at the both ends of 
the chain leading to the 4-fold quasi-degeneracy as far as 
j < 1. In this case, the energy gap G{N,j,S) decreases 
exponentially with N. On the other hand, for large (5, 
there remain no residual spins. Therefore the energy 
gap G{N, j, (5) remains finite in the thermodynamic limit. 
Thus the product NG{N, j, S) decreases (increases) with 
N for 6 < 6c {S > Sc)- Thanks to this situation, we can 
accurately determine the critical point using the phe- 
nomenological renormalization group method. 

If we search for the critical point varying S with fixed 
J, the phenomcnological renormalization group equation 
for the finite size critical point Sc{Ni,N2) reads,!!!! 

NiG{Ni,j,Sc{Ni,N2)) = N2G{N2,j,5c{Ni,N2)). 

(3.1) 

We also have the corresponding equation for jc{Ni, N2) 
if the roles of j and S are interchanged. The finite 
size critical points are obtained as the intersection of 
NiG{Ni,j,S) and N2G{N2,j,S) as shown in Fig. 2. 
Three intersections for (iVi,iV2) = (12, 16), (16, 20) and 
(20,24) are represented by double circles. 

These values are extrapolated to A'^ ^ c» assuming 
that the finite size correction to the critical value is pro- 
portional to ( ^1+^2 ^-i/i^ where v is the critical expo- 
nent of the correlation length or, equivalently, the energy 
gap. This assumption is based on the finite size scaling 
hypothesis.EB' At j = 1 and (5 = 0, the exponent v is 
known to be 2/3 with the logarithmic corrections to the 



pure power-law behavior H' a Assuming that i> remains 
constant over the whole phase boundary, we take v = 2/3 
in the extrapolation of the critical point. This assump- 
tion is plausible because the system remains SU{2) sym- 
metric over the whole phase boundary. 

To check this assumption numerically, we estimate the 
value of the exponent v from the numerical data using 
the following formula of Spronken et alB> , 
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(3.2) 



which takes into account the logarithmic corrections 
due to the SU{2) symmetry of the isotropic Heisen- 
berg model. Here G'{N) denotes the derivative of 
G{N,j,S) with respect to j or S at j — jc{Ni,N2) or 
S^SciNi,N2). 

Fig. 3 and Fig. 4 show the examples of the extrapola- 
tion procedure of the critical points and critical exponent 

for (5 = 0.5. The cross is the critical point jc — 0.889 
and exponent v ^ 0.670 in the thermodynamic limit. 
The ^-dependence of v is shown in Fig. 5. The horizon- 
tal line is I' = 2/3. It is verified that v ~ 2/3 over the 
whole phase boundary except for the neighbourhood of 
J = and S = 1 where the numerical accuracy becomes 
worse. In this region, however, the phase boundary and 
the critical exponent can be determined by the pertur- 
bation theory with respect to j as follows. 

For small j, the spins connected by the 1 -I- (5-bonds 
form strong singlet pairs. Between the spins S41+1 and 
84,1+4, the effective coupling jg-g- is generated mediated 
by the quantum fiuctuation in the singlet pairs on the 
1 + (5-bonds. This is calculated as, 
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1.-2 

2-J 

l+S 



(3.3) 



up to the second order in j. When 1 — (5 = jgg, the effec- 
tive Hamiltonian is the uniform spin-1/2 AFHC which 
has a gapless ground state. Thus the phase boundary is 
given by, 



2' 



(3.4) 



and the gap exponent v — 2/3. 

The ground state phase diagram of the present model 
is shown in Fig. 6 by the solid line. The perturbation 
result is shown by the dashed line. The numerical result 
is consistent with the perturbation theory. 

On the other hand, TakanoO^ has derived the gapless 
line of the present model as S = (1 — j)^^^ using the 
mapping onto the nonlinear (7-model. This is in qualita- 
tive agreement with our result as shown in Fig. 6 by the 
dash-dotted line. 

Each region of the phase diagram is illustrated as fol- 
lows. When j is positive and close to unity, the phases 
inside and outside the solid line are the S = 1/2 dimer 
phases with different parity. The phase boundary is 
the gapless line where the energy gap disappears due 
to the competition of dimerization and quadrumeriza- 
tion. When j is negative, the ground state is the Hal- 
dane state inside the solid line and 5 = 1 dimer state 
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outside the solid line. The 5=1/2 dimer phase is thus 
connected with the Haldane phase continuously inside 
the solid line even in the presence of quadrumerization. 
On the other hand, the connection between the opposite 
parity S = \/2 dimer phase and the 5 = 1 dimer phase 
outside the solid line is interrupted by the point 5 — \ 
and i — Q. At this point the ground state consists of 
isolated dimers on the 1 + (5-bonds and free spins on the 
4Z-th and 4/ + 1-th sites. In the presence of small but 
finite j, however, the spin structure does not depend on 
the sign of j because the effectiv e co upling j^^^ is propor- 
tional to p according to Eq. (3^). Therefore we may 
conclude that the 5 = 1/2 dimer phase with j > 1 and 
the 5 = 1 dimer phase also belong to a single phase. 

§4. Analytic Results for A = 

For A = (ATF- model), our model reduces to the 
half-filled noninteracting spinless fermion system by the 
Jordan- Wigner transformation. After the Fourier trans- 
formation, the single particle excitation spectrum is de- 
termined by the eigenvalue of the 4x4 matrix. 
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(4.1) 



where k is the momentum of the excitation. The eigen- 
value is given by. 



where 



£±(fc) = 



£ = ±e±(fc). 



{f + [l + 5')] 



(4.2) 



± V p cos2 2k + j2<52 sin^ 2k + 



1/2 



(4.3) 



In the half filled case, negative branch is filled and the 
excitation energy is determined by £±{k). The energy 
gap G is determined from the minimum of £±{k) as. 



G 
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Setting G = 0, the critical line can be exactly calcu- 
lated as. 



f + 6' = l, 



(4.5) 



and the critical exponent = 1. The phase boundary is 
a circle as shown in Fig. 6 by the dotted line. 

When j is positive, the phase diagrams of the isotropic 
case and that of the XY case look similar. For nega- 
tive values of j, however, these two models behave quite 
differently. In the isotropic case, the Haldane-like phase 
extends to j — > — oo for small S while it end up at j = — 1 
in the XY case. The phase at j — > — oo is continuously 
connected with the 5=1 dimer phase. This can be 



understood in the following way. 

The coupling on the j-bond can be rewritten as, 



jS2l-lS2l 
jS2l-lS2l 



CV qV 



-J(A- 



'^^21-1^21) 

1 



If we regard S21-1 + as a single 5=1 spin operator, 
the last term corresponds to the single site anisotropy 
term which has large positive value for j — > — 00 and 
< A < 1. Therefore the ground state at j — > —00 is 
the large-I? like phase in the XY case. Considering that 
no phase boundary exists betweeiilhe large- 13 phase and 
the dimer phase in 5 = 1 AFHCjiS-* it is reasonable that 
we find no phase transition as a function of S for large 
enough negative j in the present model. 

§5. Summary and Dicussion 

The ground state phase diagram of the dimerized and 
quadrumerized spin-1/2 Heisenberg chain is calculated 
by numerical diagonalization method. In the isotropic 
case, the critical points are determined using the phe- 
nomenological renormalization group and the finite size 
scaling hypothesis. The gap exponent v is estimated to 
be close to 2/3 over the whole phase boundary taking 
the logarithmic corrections into account. It is suggested 
that the phase transition of the present model belongs to 
the same universality class as the 5=1/2 dimerization 
transition as expected from the symmetry of the system. 
The numerical results are consistent with the perturba- 
tion theory for 5 ~ 1 and j ~ 0. 

It is found that the Haldane phase is connected with 
the 5 = 1/2 dimer phase even in the presence of 
quadrumerization S. The 5=1 dimer phase is found 
to be connected with the 5 = 1/2 dimer phase with op- 
posite parity. At first glance, the spin structures of the 
last two states might appear different. The difference 
of these two spin structures is, however, local and can 
be transformed via a local reconstruction of 4-spin sin- 
glet without closing the energy gap. Therefore these two 
states are continuously connected with each other. Ac- 
tually, no evidence of the phase transition is observed in 
the numerical results. 

In the isotropic case, our problem has always the 
SU{2) symmetry and this fixes the value of the criti- 
cal exponents. If we violate the SU (2) symmetry by in- 
troducing the anisotropic exchange interaction, we may 
expect wider variety of critical phenomena. Consider- 
ing the big difference between the phase diagrams of the 
XY case and the isotropic case, it must be of interest 
to investigate the intermediate region < A < 1 in de- 
tail. The study of this problem is in progress and will be 
reported elsewhere. 

Another way of breaking the SU (2) symmetry is to 
introduce the magnetic field. The effect of the magnetic 
field is of special interest related with the magnetiza- 

£ffarded as the field induced 
l3 This is also under inves- 
tigation and will be reported elsewhere. 
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Fig. 1. The schematic spin structure for (a) j < I, 5 < Sc, (b) 
j < 1, S > 6c and (c) S = 0, j > 1 



Fig. 2. The j-dependence of NG(N,j,S) with S = 0.5 for 4Af = 
12, 16, 20 and 24. The intersections (double circles) are the finite 
size critical points. 



Fig. 3. The extrapolation procedure of finite size critical point jc 
for 5 = 0.5 



Fig. 4. The extrapolation procedure of finite size critical expo- 
nent u for 5 = 0.5 



Fig. 5. The <5 dependence of the numerically obtained critical ex- 
ponent u. Filled (open) circles represent the values for j > 0{j < 
0). 



Fig. 6. The phase diagram of the isotropic model (solid line) and 
the XY model (dotted line). The dashed line is the result of 
the perturbation calculation for 5 ~ 1 and j ~ 0. The dash- 
dotted lips is the phase boundary obtained from the nonlinear 
a modeLEJ 
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